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Abstract 

We consider the magnetic properties of high T c cuprates from a gauge the- 
ory point of view, with emphasis on the underdoped regime. Underdoped 
cuprates possess certain antiferromagnetic correlations, as evidenced, for ex- 
ample, by different temperature dependence of the Cu and O site NMR re- 
laxation rates, that are not captured well by slave boson mean field theories 
of the t-J model. We show that the inclusion of gauge fluctuations will rem- 
edy the deficiencies of the mean field theories. As a concrete illustration of 
the gauge-fluctuation restoration of the antiferromangetic correlation and the 
feasibility of the 1/N perturbation theory, the Heisenberg spin chain is an- 
alyzed in terms of a 1+1D U(l) gauge theory with massless Dirac fermions. 
The 1/iV-perturbative treatment of the same gauge theory in 2+1D (which 
can be motivated from the mean field 7r-flux phase of the Heisenberg model) 
leads to a dynamical mass generation corresponding to an antiferromagnetic 
ordering. On the other hand, it is argued that in a similar gauge theory with 
an additional coupling to a Bose (holon) field, symmetry breaking does not 
occur, but antiferromagnetic correlations are enhanced, which is the situation 
in the underdoped cuprates. 

PACS numbers: 74.25.Ha, 74.72.-h, 75.10. 
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I. INTRODUCTION 



The essence of the physics of high T c cuprates boils down to the problem of how to treat 
the dual nature of the electrons which form local moments in the insulator, and yet make 
up a Fermi surface when doped with ~ 15% holes. This problem was brought into sharp 
focus with the discovery of anomalous properties in the physics of the "underdoped" region 
which lies between the antiferromagnetic insulator and the optimally doped superconduc- 
tor. How does the Fermi surface evolve from small hole pockets near k = (±7r/2, ±vt/2) in a 
slightly doped antiferromagnet to the full Fermi surface obeying Luttinger's theorem in the 
optimally doped materials? What are the magnetic properties in this intermediate doping 
region? Experimentalists have already answered a substantial part of these questions. In 
particular, the angle-resolved photoemission spectroscopy (ARPES) has shown the existence 
in the normal state of a gap with the same anisotropy as the d-wave gap of the superconduct- 
ing state [JTJ] . Low-lying excitations are observed along a patch near (±7r/2, ±7r/2) ["Fermi 
surface segments"], but the Fermi surface, in the veritable sense of the word, does not exist. 
The ARPES results might have been accepted without much grudge simply as a plausible 
interpolation between the antiferromagnet and the optimally doped superconductor, had 
our understanding of metals not been so entrenched in the Fermi liquid theory; the notion 
of a metal without a Fermi surface is a serious embarrassment. At the same time, gap- 
like suppression of spin excitations are seen in NMR 0: the Knight shift and spin-lattice 
relaxation rates all decrease with decreasing temperature below certain temperatures. 

Gaplike features in the underdoped cuprates might remind us of the spin liquids — the 
liquid of spin singlets. Yet the devil is in the details, and the underdoped cuprates deviate 
significantly from "gapped spin liquids" like spin 1/2 ladders and integer spin chains. The 
latter have an even number of spins per unit cell and the ground state is a spin singlet. They 
are characterized by a clear gap to the lowest triplet excitation that is inversely related to 
the correlation length; this gap can be seen in inelastic neutron scattering. The magnetic 
responses like uniform susceptibility and the NMR relaxation rate as a function of tempera- 
ture have activated behaviors. Many of them can be satisfactorily described in terms of the 
"quantum disordered" phase of the nonlinear sigma model C = ^(d^n) 2 ,g > g c || which 

can be also understood in terms of the CP N ~ X model involving bosonic "spinons" (z fields) 
- these spin 1/2 excitations are said to be confined, as they do not appear in the basic 
physical spectra In other words, a description of the excitations in terms of fluctuating 
spin 1 objects is most natural for them. 

On the other hand, in the underdoped cuprates, inelastic neutron scattering does not 
show a "full gap" ||. The decrease of the Knight shift with decreasing temperature looks 
more like a power-law. On the whole, the magnetic excitation spectrum of the cuprates 
seems to display a curious mixture of singlet and antiferromagnetic correlations. There are 
evidences for antiferromagnetic correlations from the Q-space scan of neutron scattering 
cross section, and from the difference in the temperature dependence of the NMR relaxation 
rates : the Oxygen \jT\T (which has little contribution from spin excitations near wave 
vector Q = (vr,7r)) monotonically decreases with decreasing temperature, while the Copper 
\jT\T (which weighs (tc,tt) spin excitation strongly) increases with decreasing temperature 
until around 150K and then falls 0. 

The ant if erromaget- singlet debate has enormous ramifications for theories of high T c 
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cuprates (for a succinct review, see Ref. ||). Some theorists have advocated the picture of 
Fermi liquid quasiparticles exchanging "antiparamagnons" for understanding the anomalous 
normal state properties and the superconductive pairing. Such a view is not shared here. 
Instead of viewing the antiferromagnetic fluctuations as the cause of superconductivity in a 
BCS-like scenario, in this paper we would rather regard them as a residual but important 
consequence of local repulsive interactions that lead to superconductivity in the presence of 
doped holes, a part of Nature's conspiracy to find a compromise between a magnetic ground 
state and an itinerant metallic state. This line of thinking goes back to Anderson's seminal 
1987 paper || on the resonating valence bond (RVB) theory, in which he reasoned that 
doped holes may propagate coherently in a liquid of spin singlets. 

Theoretical attempts to realize Anderson's RVB picture are based on strong coupling 
models, such as the one band Hubbard model or the t-J model ||. These models are 
considered to contain some essential physics of the cuprates at appropriate parameter values 
U/t or J/t. The t-J model, the simpler of the two, captures in a transparent way what is 
believed to be the basic physics, namely the competition between the magnetic exchange 
and the derealization energy of holes. The no-double-occupancy constraint in the t-J model 
can be taken care of by writing the electron operator as a composite of a neutral fermion 
(spinon) and a spinless boson (holon) [c] CT = f}j)i\ and demanding each site be occupied 
by either a fermion or a boson (J2 a AvAr + b\bi = 1). The theory then contains four- 
particle interactions which can be decoupled by introducing "mean fields" like Xij — (fLfjv)> 
\j — (At At ~ At /it)> an d Vij = (bjbj). Within the mean field approach, Kotliar and Liu, 



and Fukuyama and coworkers [10] have studied the phase diagram of the t-J model. At 



low doping (and below some temperature scale), it was found that the phases in which the 
fermions are paired into d-wave singlets (A^+a, = —A iti+y ^ 0)are favored. Depending on 
whether the bosons are condensed, they could be superconducting (SC phase) or normal 
("d-wave RVB phase"). 

As noted by Rice [|nj and others, the fermionic mean field theory captures some impor- 
tant features of the spin gap phenomena in the underdoped cuprates that refuse clear-cut 
characterization in terms of a well-defined correlation length and relaxation times. The 
theory describes some kind of quantum spin liquid, but unlike gapped spin liquids, there 
is a particle-hole (spinon-antispinon) continuum, which would create some spectral weight 
for magnetic excitations at arbitrarily low energy. More specifically, the Dirac spectrum 
(e(k) = v |k|) of the fermionic quasiparticles in the d-wave RVB phase gives the Knight shift 
K ~ T and the Oxygen site NMR relaxation rate 1/T\ ~ T 3 , in rough agreement with 
experiments in the underdoped cuprates. Moreover, the absence of the gap in the charge 
response (for example, the in-plane optical conductivity) could be explained simply, since 
the spin and charge degrees of freedom are separated, i.e. the spin is carried by fermionic 
spinons while the charge is carried by bosonic holons. 

A Dirac-type spectrum as in the d-wave RVB phase was also found by Affleck and 
Marston who considered the 7r-flux phase [0] as a possible spin liquid ground state of the 
cuprates. It turned out that at half filling the d-wave phase with |Ajj| = \xij\ is equivalent 
to the flux phase, due to a local SU(2) symmetry [|13f| . Wen and one of us reasoned that 
this symmetry might still be a pretty good (and important) symmetry at small dopings, and 
came up with a slave boson theory that respects the SU(2) symmetry even away from half 
filling by introducing an SU(2) doublet of slave bosons, hoping to get a better description 
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of underdoped cuprates WM. In this theory, the mean field corresponding to the "spin 



gap" phase of the underdoped cuprates was identified as the "sFlux phase" which can be 
considered a combination of the d-wave RVB phase and the staggered flux phase |T3| of the 
U(l) theory. This phase also has fermions with a Dirac spectrum, but in contrast to the 
d-wave RVB phase of the U(l) mean field theory [|16[], the fluctuations around this mean field 
include a massless gauge field which is expected to affect strongly the magnetic and transport 
properties of the system E7| . With the inclusion of a residual attraction between bosons and 
fermions, the sFlux phase was shown to reproduce the gross features of the ARPES, such 
as the Fermi surface segments near (±7r/2, ±7r/2) and a large gap at (tt, 0). 

Despite the successes, the mean field treatments (both the SU(2) theory and its pre- 
dessesors) are unsatisfactory in several respects. For example, it is not clear how the spin 
gap phase is connected to the Neel ordered phase at zero doping. The mean field ansatz 
loses a lot of antiferromagnetic correlation; within the mean field theory, the Copper site 
1/TiT has a similar behavior as the Oxygen site 1/TiT, in disagreement with experiments 
0. Attempts to fix the problem by some kind of RPA cannot produce, in a natural manner, 
different temperature scales which mark the decrease of the Copper and Oxygen site relax- 
ation rates [Tj| . Another serious question is the role of gauge fluctuations around the mean 
field solution which had not been studied carefully so far. In fact, a strong gauge fluctuation 
might destroy the mean field picture altogether, in which case we have to re-identify the 
elementary excitations of the theory. 

In this paper, we look into these questions. The basic point is that the gauge fluctuations 
ignored at the mean field level strongly enhance antiferromagnetic correlations. The gauge 
fluctuation could be so strong that the elementary excitations of the mean field theory 
disappear completely from the low energy spectrum. This is believed to be the case with 
the "U(l) 7r-flux phase" (I2j description of the undoped cuprates (antiferromagnet), which 
is a theory of massless Dirac fermions coupled to a U(l) gauge field (massless QED3). The 
QED3 can be treated in the 1/N perturbation theory. For physical N(= 2), a dynamic 
mass generation and spontaneous symmetry breaking corresponding to Neel ordering would 
occur [19 1, while for large enough N, the theory would still describe some kind of a spin 
liquid. The true low energy excitations of the symmetry-broken case are recognized as the 
Goldstone bosons- "mesons" which are a bound state of a particle and an antiparticle (spinon 
& antispinon). 

We note that the local SU(2) symmetry at half filling can be utilized to write a theory 
of the undoped system in terms of massless Dirac fermions coupled to nonabelian (SU(2)) 
gauge fields [fL3| , |i~7[[ . Since the U(l) theory (at least initially) has one massless gauge field 
while the SU(2) theory has three, they look quite different. Nonetheless, as long as the gauge 
fluctuations are treated exactly, the two theories should lead to the same physics, namely 
the Neel ordering. 

As regards the underdoped cuprates, we believe that the most promising starting point 
is the sFlux phase of the SU(2) theory ||17|| , which has massless Dirac fermions coupled to 
a massless U(l) gauge field as in the Affleck-Marston flux phase ||12[| , but also has nonrel- 
ativistic bosons (ho Ions) coupled to the same gauge field. We shall argue that, due to this 
additional coupling, the gauge fluctuations will not destroy the essential validity of the mean 
field picture, though the picture of antiferromagnetic spin excitations will be modified (im- 
proved). The coupling to the bosons would result in the screening of the time component of 
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the gauge field which will prevent Neel ordering. The gauge field will nevertheless mediate 
an attraction between spinons and antispinons and try to create a bound state with mo- 
mentum ~ (tt, tt), but due to the particle- hole continuum, this will appear only as a broad 
resonance. This can be viewed as a Goldstone boson precursor mode that comes down in 
energy as the transition is approached (as the boson density is reduced). The recent neutron 
scattering in underdoped cuprates which sees a broad peak in Q « (tt, tt) magnetic response 
whose energy scale is roughly proportional to doping might be consistent with this point of 
view [22]. We shall also discuss the issue of confinement, as there are lingering questions 
about the fate of "spinons" in the case of strong coupling gauge theories. 

In order to illustrate some aspects of the foregoing ideas more concretely (in particular 
the gauge-fluctuation restoration of antiferromagnetic correlation and the feasibility of a 
U(l) gauge theory description of quantum antiferromagnets) , we'll first reexamine the well 
known spin half chain from the point of view of the Schwinger model. 



II. LESSONS FROM SPIN CHAIN 

We begin with a discussion of the Id Heisenberg model as an example where the idea 
of mean field theory plus gauge fluctuation can be applied to an exactly soluble model. A 
number of authors have noted that what emerges is the Schwinger model with two flavors 
25|, and abelian p4j and nonabelian [2l| bosonization methods have been applied to solve 



for the continuum limit of this model. Here we solve the model using a slightly different 
technique, i.e., chiral rotation, which focusses on the role of gauge fluctuations as a way of 
correcting the underestimate of the antiferromagnetic correlation in mean field theory. We 
then show that qualitatively similar results are achieved in 1/N perturbation theory of the 
gauge fluctuations. This gives us encouragement that a similar approach may be reasonable 
in the two-dimensional case. 



A. RVB theory of spin one-half chain 

The Heisenberg model (H = JJ2<ij> Sj • Sj) in the fermion representation of the spin 
can be written 

H = -- fjafiafipfjP W 

Z <ij> 

with the constraint Y^aflafia = 1- The 4- fermion interactions and the constraint can be 
handled by the introduction of a Hubbard-Stratonovich field and a Lagrange multiplier, 
which gives 

H = -J/2 £ (XijfLfja + h.c.) + • £ Xiiflfia - 1). (2) 

<ij> i 

Within the mean field theory, Xij — Xi ^ = 0> hence the mean field hamitonian in the k-space 
is 
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H mf = -x-J ]T cos(k)flj ak (3) 

k 

= -XJ^k COs(k) (flj ak - flk-Ja,k-ir) (4) 
= -X-J^'k COs(k)(fl e J aok + ft ok faek)- (5) 

Here H' k denotes sum over the magnetic BZ, say < k < ir, and f ek , fok are even and odd site 

operators (f ek = ^(f k +fk-*) = /Steven e ikj fj, fok = ^{fk~fk-n) = £jodd e ikj fj). 
Linearizing around k = n/2 + k! , we arrive at the continuum hamiltonian 

H = -Jdk'^ a (k')cx l k'Mk'), (6) 

where ip a = ^ ^ and o\ = ^ ^ J ^ is a Pauli matrix. This is just the hamiltonian of free 

Dirac fermions (we have set the velocity of the fermions =1). The corresponding (Euclidean 
space) lagrangian is 

L = ^al^d^a, (7) 

where ip = V^7o> an d fx — 0,1, and the 7 matrices are 

7o = 0-3, 7i = -°2- (8) 
In 1+1D, we define 75 matrix as 75 = — ?7o7i = o\, which has the property 

{75,7^}=°^ VTi/ = »7m75, (9) 

[e^y is the antisymmetric tensor with eoi = 1). 

Including the fluctuations around the mean field (the fluctuations of \ and the phase 
of Xij) amounts to coupling the fermions to a U(l) gauge field by the minimal prescription. 
Hence the continuum version of {the mean field + fluctuations} is the two flavor Schwinger 
model 

L = i>al^(d^ - ia^ipa. (10) 

The apparent gauge coupling (bare coupling) is infinitely strong, as there is no kinetic term 
for the gauge field. 

Above lagrangian is obviously invariant under the global SU(2) transform (spin rotation 
symmetry) 

lp a -> (exp(^V)) Q /#/3 (11) 

where r l , I = 1, 2, 3 are Pauli matrices (belonging to a space different from that of 07). The 
lagrangian is also invariant under the "chiral transformation" 

ip a -> exp(-?6»7 5 )^ a . (12) 

This "chiral symmetry" is explicitly broken by higher derivative terms ignored in taking the 
continuum limit, like 
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(13) 

fk— ki where < k < ir, 



If we label the states by the left and right movers, fn~fk, fz 
the chiral transformation corresponds to fn — > jRt~ %e , fi — > fh^ 6 - 

We now consider spin correlation functions at the mean field level (i.e. ignoring gauge 
fields). The spin operators in the continuum has two contributions (uniform & staggered): 

To evaluate the spin correlation function 

= (4>uoMx)4>2ioMQ)) + (-if 1 (4>iMx)4>2M0)) 

(t = (r 1 ±ir 2 )/2), we need the fermion Green's function 

G a/3 {x) = (^a(aj)^(O)) = G(x)5 a p 
which can be obtained from the momentum space Green function G(k) = —ik^/k 2 as 



G(x) = 7^ 



d 2 k e 



—ik-x 



d 



d Xfl J (2vr) 2 k 2 



2nx 2 



(14) 

(15) 
(16) 

(17) 



Here and from now on, unless otherwise specified, we use the usual field theory notation: 



k = (k ,ki),x = (xq,xx) [italics denote space time vectors]; ary 
Using Wick's theorem, we have 



XfjTfn] x 



Xq ~\~ x ^ , etc. 



(S + (x)S~{0)) = -tr 7 , T [G(s)7or + G(-x)7 r-] - (-l) Xl tr 7>T [G{x)r + G{-x)T- 



2tt 2 
1 

4^2 



Xq X 



(Xq + x 2 ) 2 



+ 



^y>2 I ^y>2 

Xq "T X\ 



+ 



+ (-1 



,x 1 



x_ 



X—X-\- 



(18) 



(x± = Xo ± ixi; tr 7)T denotes trace over both the 7 and r spaces (spinor and spin spaces)), 
which does (and should) equal the (S z (x)S z (0)) correlation function in the XY model |26| . 
The equal time correlation function (S(xi) • S(0)) behaves as 

( S (^)- S (°)) = ^4 1 ? (( " ir " 1) ' (19) 

(the spins on the same sublattice are not correlated at all, while the correlation among spins 
on different sublattices are decaying algebraically as 1/x 2 ). This peculiar behavior, which 
was derived by Arovas and Auerbach in the lattice version of the fermionic mean field theory 
P7J and agrees with Bulaevskii's Hartree-Fock treatment of the Jordan- Wigner fermionized 
Heisenberg model ||28|| , is viewed as a pathology of the mean field theory: we have lost a 
substantial amount of antiferromagnetic correlation. 
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B. Schwinger model 



We now consider the effect of gauge fluctuations. It is natural to expect that the in- 
clusion of gauge fluctuations will improve the mean field picture. The time component 
of the gauge field can be regarded to originate from the Lagrange multiplier field (for no- 
double-occupancy); this corresponds to Gutzwiller-projected (half-filled tight binding) Fermi 
surface, which is known to be a pretty good description of Id antiferromagnet [[2^,|3(]] • 

As mentioned earlier, our theory with fluctuations is a Schwinger model. For reasons 
that will become clear shortly, we consider a slightly more general case of iV-flavors: 

DipDipDa^ exp(— S), 

N _ j 

S = d 2 xJ2 ^al^dp - ia^)ip a + -t~1 F 1^ {<? = oo). (20) 

J a=l 4e 

The physical case is iV = 2; general (even) N corresponds to an SU(iV) antiferromagnet. 
Integrating out the fermions gives 

Z = J Da^ exp(iV Tr ln(l - iQ^aJ) (21) 

where Q(x, x') = ( r y fl d^ 1 5(x — x'), and Tr denotes traces over the spinor space and the 
position space (Tr=tr/ d 2 xd 2 x' '...). The logarithm can be expanded, giving 

Z = J Da^exp (~ J d 2 xd 2 x'a fl (x)Il^(x - x')a u (x')^j , (22) 

where U^x) = — N tr[G(x) , y lJi G(— x)j u ]. Note that the beyond- Gaussian terms (like 



r fj.upSO'iJ.O'vO'pO's) are all zero; the proof can be found, for example, in Ref. |31 . 

The polarization function H^^q) (in the momentum space) contains a divergence that 
has to be regulated using gauge invariant schemes, like the dimensional regularization or the 
Pauli-Villars regularization. Relegating the details to Appendix A, we have 

IW*) = - (V - , (23) 



7T 



k 2 



which means that the gauge boson acquires an infinite mass (= e^N/n). The transversality 
of Eq.|^ guarantees the conservation of the current j M = Tpal^a'- 

qMq) = Q l i(^u(q)a v (q)) = -> <9 M j M = 0. (24) 

On the other hand the current j5 At (= i^pal^l^a = —t^uju) associated with the chiral sym- 
metry (Eq.|i~2"D is not conserved: 

iN . iN 

Qfijb/j, = Qfi^uvi^-upO'p = ^auQri^i/ * dujz>u = — £nvFn U . (25) 

7T Z7T 

This result, the so-called axial anomaly, can be regarded as either a consequence of or a 
condition for gauge invariance. Equation (25) is in fact familiar to the solid state physicists 



32|| . In an electric field, the equation of motion for the crystal momentum state is 
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(26) 



In one dimension an electric field causes a shift of the occupation between left moving and 
right moving states. Thus the left and right movers are in fact connected and their densities 
are not separately conserved. It is easy to see that Eq.(26) is consistent with Eq.(25). 

The exact spin correlation functions can be evaluated with the use of "chiral rotation" 
j33H36|. In this approach, the gauge field is written as the sum of a div-free part and a 
curl-free part: 

a M = V^A + (27) 

The transform 



ip a -> ip' a = exp(-7 5 6 l a - i9 b )i> a (28a) 
4 -> = ?Pa exp(- 75 # a + iO b ) (28b) 

decouples the gauge field from the ip' fermions. We note that corresponding to Eq.(28b), we 
have 

$i->Vtfexp(750a + i0&) (29) 

because 7 and 75 anticommute. Thus ip^ipa = ip'd^'a an d the transformation is unitary. 
Here ip a and ip^ are treated as independent Grassmannian variables and exp(7 5 # a ) should 
not be thought of as an amplitude transformation. Alternatively, we point out that in real 
time (as opposed to imaginary time used here), the transformation to cancel the a M field will 
indeed be a phase rotation and takes the form ip a — > ip' a exp(—i'y 5 8 a + i6 b ). That this must 
be the case is evident from Eq.(27), because ao = d\6 a changes from real to imaginary upon 
going from Minkowski to Euclidean space. Unlike the usual gauge transformation, zoq/l is 
cancelled by difi and iaifi is cancelled by dofi under the chiral rotation, so that 9 a must 
also change from real to imaginary. 

Even though the chiral rotation is unitary, it does not leave the Grassmann measure 
invariant, due to the fact that fa and fi are not truly independent, as noted earlier. The 
jacobian exp(J) for the change of measure 

DipDip = Dij)'D$ exp(J) (30) 



can be found easily using the axial anomaly condition Eq.|25|(see Appendix B for details): 



N r 

J=~ d 2 x(d,e a ) 2 - (31) 

Z7T J 

Alternative derivation can be found in Ref. |36 ]. In terms of the new fields, the functional 
integral is 

Z = J DftDiP'D9 a exp - jd 2 xU' a lM' a + |^aA) 2 ) . (32) 
These are "free" fields, and now the spin correlation function (Eq.|l5D can be evaluated easily: 
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(^i7o^2 0)^270^1 (0)) = (^'170^2(^)^270^1 (0)) 

1 

= 2tt 2 (xl + xj) 2 ' 
(^2(^)^2^1(0)} = (^e 2 ^^(^e 275 ^0)^ 

_ 1 c 2((e a (x)~8 a (0)) 2 ) 



2tt 2 x 2 

x 2 (x 2 y-v N ' 



(33) 



where we have used 



((^W-^(0))^-|/ A |^i(^-i) 

= ^M*W), (34) 

(A is a UV cutoff originating from the lattice theory, and C is a nonuniversal constant that 
depends on high energy details A). In the physical case (N — 2), we then have 

(S(*)-S(0))~(-ir-=L= (35) 

which agrees with the more accurate result [37]] up to a ln : ^ 2 (x 2 ) factor. The log factor, 
not captured by the Schwinger model, must be due to terms ignored in our derivation from 
the lattice theory, e.g., the amplitude fluctuation of the RVB field; this is analogous to 
the bosonization theory of Heisenberg model, in which the Umklapp processes give rise to 
logarithmic correction (prefactor) to the power law ||38|| . The simple correlation function of 
Eq.|55| is known to capture the low energy (temperature) properties of the Heisenberg spin 



chains quite well [39 



C. additional remarks 

Before moving on to the perturbative treatment of the same theory, we briefly compare 
the gauge theory approach with other approaches. Perhaps the best known treatment of 
the Heisenberg spin chain relies on the Jordan- Wigner transformation (fi = e^S 1 " , (pi = 
7rX)}=i St^r)' However, unlike the X-Y model, this does not lead to a theory of free 
fermions, but to a theory with a 4-fermion interaction with coupling constant of order unity, 
which is then treated by bosonization methods |38| . A drawback of the Jordan- Wigner 
approach is that the SU(2) spin symmetry is easily lost, and the correct exponent (= —1) 
of the staggered spin correlation function has to be determined rather indirectly. 



Haldane [40[ , on the other hand, proposed to analyze the Hubbard model, exploiting the 



"equivalence" between the spin sector of the Id Hubbard model 

H = -\{JU) 1 ' 2 £(ie i+v + h.c.) + ^(n it - i)(n a - \) (36) 
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and the Id Heisenberg model, that holds even in the weak coupling limit. The Hubbard 



model can be treated by a bosonization method [40,26 that respects the SU(2) symmetry at 
all stages, and the correct exponent of the spin correlation function can be obtained directly. 
We can easily adapt Haldane's approach to the Schwinger model, which as we point out, 
is equivalent to the RVB mean field theory plus gauge fluctuations. This line of approach 
was recently given in ref. [24]. Instead of treating the gauge fluctuation by chiral rotation, 
we first integrate out the gauge field in Eq.(20), which simply enforces the constraint of no 
double occupation, i.e., the charge fluctuation is completely suppressed. Now we can treat 
the fermion problem by bosonization in the standard way. We write 

(fsR(Xl)JsL(Xl)) = e i ^Mxi)^ e iy^e a (x 1 )+ik F x 1 ^ e -i^e a { Xl )-ik F x^ ( 37 ) 

for s = spin up or down, and T = ^= (9 P + 9 a ), Q x = ^ (9 P - 9 a ), fa = ^= (<f> p + <f) a ), fa = 

(0 P — fa), where d Xl 9 p ^ is proportional to the charge (spin) density and is the conjugate 
variable to <f)p( a ). In the effective Lagrangian, spin and charge fluctuations separate. For free 
fermions 9 a and 9 P contribute equally to the 2fcp spin-spin correlation functions, yielding 
Eq.(19). In our case the charge fluctuation is suppressed and we may ignore the 9 P degrees 
of freedom. The 2k f spin correlation function decays with an exponent exactly one half than 
that of the free fermion, and we recover Eq.(|33]). This line of reasoning has a slight advantage 
over Haldane's original treatment in that the energy scale J for the fermion bandwidth is 
correctly produced. 



D. perturbation theory 

We now examine the gauge theory in terms of the perturbation theory in 1/N. The 
key point is that the nature of the perturbative correction to the mean field results is 
very different for the uniform part and the antiferromagnetic (staggered) part of the spin 
correlation: While Q = n response is strongly affected by perturbative correction, the Q = 
response receives no correction at all (the absence of any correction is a special feature of 
the Id). 

The leading 1/N correction to the spin correlation functions can be straightforwardly 
evaluated. They are represented by the Feynman diagrams in Figjl] and Fig.^|. Within the 
usual Faddeev- Popov scheme of gauge fixing (the introduction of the ^r(d ■ a) 2 term to the 
lagrangian), the gauge propagator (represented by the wiggly line) is given by 

D„, = M = £ (V ~ S f)+ • (38) 

We choose the Landau gauge (A = 0) which is a natural choice, since in this gauge no infrared 
divergence occurs in the perturbation theory. The fermion propagator G(p) = {ip'y)" 1 is 
represented by a solid line. The fermion-gauge vertex is simply Z7 M ; the external current 
vertex is 70 for the uniform part (represented by a square) and 1 for the staggered part 
(represented by a circle). Of course a trace is taken over the fermion loops. 
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FIGURES 




a) b) c) 

FIG. 2. Leading 1/N correction to staggered spin correlation. 
In 1+1D, the transverse projector has a special property 

which, together with the "Ward identity" 

G(p + q) qi G{p) = i{G{p) - G(p + q)), (40) 

simplifies the algebra substantially (Note •y^e^pPp = —^775). For example, the diagram |I]b 
is 

1 " "' ' tr[G(p + q)q' ll5 G(p + q + q , ) l0 G(p + q')q' ll5 G(p) l0 ]/q' 2 



NJ (27r) 2 (27r) 2 

= ~ / ^^{( G (P + q)-G(p + q + q')ho(G(p + q') - G(p)) lo \/q>* . 

(41) 

It's straightforward to show that sum of the diagrams |l|a+c is the same as Eq.|41], except for 
a minus sign. Therefore in the uniform channel, the vertex correction and the self energy 
correction cancel. Similar cancellation is expected at all orders of perturbation theory; the 
nonrenormalization of uniform part of the spin correlation function is quite natural, since in 
our theory 

(^170^2^270^1) oc (joJo) = n 00 , (42) 
and Eq.|22] is an exact result. 
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On the other hand, the diagrams in the staggered channel do not cancel. The sum of 
the diagrams ||a and ||c are equal to ||b, which is given by 

I / (2^2 ^ tr K g (P + Q)-G(p + q + q'))l(G(p + q>) - G(p))l]/q' 2 

Therefore, in the coordinate space, the 1/N correction is 

*P(*)G<-*)]£/ |^(e<- - 1) = J^H* 2 )- («) 

Similar (but a lot more tedious) calculation would show that 1/N 2 correction is given by 
4^2^23,2 ln(x 2 ). In other words, the perturbation series exponentiates: 

$W*(aO&lMO)> = 2^(1 + ln(a; 2 ) + 1 -{1/N) 2 \r?{x 2 ) + ...) oc ^ 2) ]_ 1/N , (45) 

giving the same result obtained from the chiral rotation approach. 



III. TWO-DIMENSIONAL UNDOPED CUPRATES 



The success of 1+1D gauge theory with Dirac fermions in describing the Heisenberg spin 
chain tempts us that a similar theory of massless Dirac fermions strongly coupled to a U(l) 
gauge field might describe a 2d quantum antiferromagnet. In fact, it is known from lattice 
gauge theories [42] that this is indeed so. In connection with high T c cuprates, Marston ||21|| . 
Laughlin |[20[| , and others have noted that the "chiral symmetry breaking" in 2+1D U(l) 



gauge theory, discussed by Pisarski []43[ and Appelquist et al. |44| in particle physics context, 
corresponds to the Neel ordering. In this section, we shall discuss this matter carefully, and 
clarify some issues related to the pattern of symmetry breaking and Goldstone bosons (the 
particle physics literature [43,44] envisions symmetry breaking pattern U(4)— ► U(2) ® U(2), 
hence 16—2-4= 8 Goldstone bosons, while the symmetry breaking pattern for Neel ordering is 
SU(2)^ U(l) which gives 3 — 1 = 2 Goldstone bosons). SU(2) gauge theories with massless 
Dirac fermions may also describe the quantum antiferromagnet [15,41)1, but we shall not 



consider this possibility because of the greater complexity of the nonabelian gauge theories. 



A. Dirac fermions and the 2d Heisenberg antiferromagnet 

A 2+ ID theory of Dirac fermions 

L = TpadftJulpa, (46) 

(fi— 0,1,2) contains fermions whose density of states behaves as ~ |e| (in a general D = 
d+1 dimensions, the density of states will be ~ [e^ -1 ). In condensed matter context, the 
paramagnetism of such fermions will result in the uniform susceptibility behaving as T d_1 . 
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This seems to have little in common with the 2d antiferromagnet, but we shall see that 
a gauge field coupled to fermions can produce the correspondence with the physics of 2d 
antiferromagnet. Such a picture can be motivated from the "7r-flux" phase mean field ansatz 
of the Heisenberg hamiltonian. The ansatz is so named since the phase of the product of 
the mean field parameter Xij around each plaquette (Imln(xi 2 X23X34X4i)) is 7r - Despite the 
"flux," the mean field does not break the parity and time reversal symmetry since the flux 
of 7i is equal to — ix. This phase has a lower energy than the BZA phase (with a large Fermi 
surface) it also has the fermion spectrum 



e(k) = v\J cos 2 {k x ) + cos 2 {k y ) 



(47) 



that roughly captures the high energy features of the undoped cuprates and the dispersion 
of a single hole in the antiferromagnet p8| , ^9| . The low energy fermionic excitations of the 
flux phase reside near two "Fermi points", ki^ = (tt/2, ±7r/2). Linearizing around these 



points gives the continuum theory [12 



where ip a i 



fale 
falo 



(i 2 



fo2o 
fa2e 



Ipaad^^lpaa, (48) 

a = 1, 2 labels the two Fermi points; e,o denote even 

we have a 



and odd sites.) Organizing the tp a i,ipa2 fields into a single spinor ip a = 
theory described by the lagrangian of Eq.E6l with 4x4 7-matrices: 



To 



03 






-0-3 



7i 



02 






-cr 2 



72 



01 
-o-i 



(49) 



Similarly to the Id case, the uniform spin 5q =0 is given by ^oTI^a^Pi while the staggered 
spin Sq_, . is given by ip a T l al3 iljp- The spin correlation function (S + (x) ■ ^"(O)) can be 
evaluated easily at the mean field level using the Green's function 



G a/3 {x) = 5 a0 G(x) 



d 3 q 



-iq-x 



dxj (2vr)3 g2 



S, 



X'J 



a/3 



4tt(x 2 ) 3 / 2 



(50) 



(as in 1+1D, the italics like x, q denote space-time vectors, and q 2 = q$ + q 2 , x 2 
We have, 



x 2 + x 2 .; 



(S + (x)S~(0)} = -tr[G(x)7oG(-x) 7 o] 



4tt 2 



• ' n 



+ 



\Xl+X2 



X" 



1 

X~ 4 



ti[G(x)G{-x)] 



The equal time correlation function is then 



(5+(0,x)5-(0,0)) 



4vr 2 x 4 



(51) 



(52) 



The mean field spin correlation function falls off as 1/x 4 , and again as in Id, the spins on the 
same sublattice are not correlated p0| , |27|| . Note, however, that due to our neglect of (vr, 0) 
contributions, Eq.[52] differs somewhat from the more accurate expressions of Refs. [ptlpTj . 
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To go beyond the mean field level, gauge fluctuation is included by the minimal coupling 
scheme, leading to the following theory: 



Z = J Dip Dip D a p exp ^- J d 3 x ^ ip a (d^ - ia^^^j (53) 

where N is a general number (in the physical case, N = 2). The absence of the j^F^ term 
means the bare coupling is infinite, but the theory is still sensible: the infrared behavior 
of the theory is well-behaved (within 1/N expansion, as we shall see), and the original 
lattice theory sets a ultraviolet cutoff. (Ultraviolet divergence can be also regulated by 
the the kinetic term -^F^ v with e 2 < oo.) As in 1+1D, the lagrangian contains more 
symmetries than the Heisenberg model. For example the theory is invariant under the 

transform ifj a — > exp(i'j4^6)il) a (where 74 = |^| ,75 = ( ^Q^;/isa2x2 unit 

matrix). Again, these symmetries are broken by higher order terms (lattice effects) which 
have been ignored. 



B. content of the gauge theory 

We now explore the physical content of the theory with massless Dirac fermions. Inte- 
grating out the fermions generates the dynamics for the gauge field (see Appendix A 
for details) 

Z = J Da^exp (~ J -^L a /1 (g)n^ I/ (g)a„(g) J , 

n^ = fv^(^-^)- ( 54 ) 

The form of indicates that, unlike the 1+1D case, the gauge field is massless, though 
the infrared behavior is not as singular as the free gauge field (i.e., L = \F*). The effect 
of the fermion-gauge field interaction can be analyzed perturbatively in 1/N. The gauge 
propagator in the Faddeev-Popov prescription is 



DM 



U_^A +X ^ (55) 



As in 1+1D (Sec. II), we choose the Landau gauge (A = 0) to avoid spurious infrared diver- 
gences. 

The fermion self-energy at the leading order in 1/N is given by 

m = 1 1 ^ 1 ^^ 1 ^ q) ~ %kl ln {%) ■ (56) 

This log divergence does not occur in the 1/N correction to the polarization function rep- 
resented by the diagrams in Fig|| These diagrams have been calculated explicitly in a 
different context by Chen et al. [pOj] . They found that the logarithmic divergences in the self 
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energy correction (Figs.[3] a+c) are cancelled by the vertex correction (Fig.[3]b); the diagrams 
sum to 



4vr 2 



(57) 



which is of the same form as the zeroth order result (Eq.^) but down by some factor 
involving 1/N. Thus, although we do not have a complete cancellation (like Id), the gauge 
field is essentially unrenormalized, except for some modification of the effective coupling 
constant (1/N — > 1/N « 1/N). In other words, Z3 (charge renormalization) « 1. Weak 
correction to the gauge propagator (despite strong self energy correction) was seen in several 
other contexts, including the half-filled Landau level (fermion Chern-Simons theory) fxTf . 
bosonic Chern Simons theories []52| , and the uRVB gauge theory |53| . This robustness is due 
to the fact that H^ v is a correlation function of conserved currents, and conserved currents 
have no anomalous dimensions (as a consequence of the Ward identiy) |54|]. The foregoing 
argument provides some optimism for perturbation theory in 1/N. 




FIG. 3. Leading 1/iV correction to vacuum poloarization. The diagrams d) and e) are zero due 
to Furry's theorem. 

Because the uniform spin correlation U u (x) = ( , 0i7oV'2(^)'027oV ; i(O)) ~ n o(^), it 
wouldn't be renormalized significantly. On the other hand, the staggered part TL s (x) of 
the spin correlation function does not involve conserved currents, and is expected to be 
strongly affected by gauge fluctuations. This difference can be more or less seen in per- 
turbation theory. The diagramatic representation of U u is the same as in the 1+1D case 
(Fig.0). The external vertices in this fermion-gauge field vertex (70). Using the 

Ward identity 



_d_ 



(58) 



we can see that the 2 overlapping divergences of Fig.[Tfb are cancelled by self-energy bubbles 
of Fig.^Ja+c. For LI S , the external vertices are 1 (unit matrix in the spinor space). In that 
case, the cancellation does not occur, and divergences develop in LI S . 

The question is, what would be ultimately the behavior of IL? Would Tl s (x) be char- 
acterized by simple power law correlations like Tl s (x) ~ l/(x 2 ) a , (a < 2) like the Id case? 
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This might be a realization of Anderson's 2d Luttinger liquid scenario [[5l|, but we feel that 
a priori there is no reason to expect so; after all, 1+1D is rather special, with all sort of 
fascinating features like the infinite dimension of the conformal group and Coleman's 
theorem |56| which prohibits spontaneous symmetry breaking. More plausibly, we would ex- 
pect a symmetry-broken phase (Neel order) or a symmetric phase with a more complicated 
magnetic correlations. 



C. spontaneous symmetry breaking 



The previous section has identified a possible antiferromagnetic instability, which corre- 
sponds to the staggered magnetization {^t 1 i/j) acquiring a definite orientation — an SU(2) 
symmetry breaking. We now examine this possibility more closely. In the symmetry-broken 
case, the fermions acquire a mass (dynamic mass generation). 

Without loss of generality, we assume that the rotation symmetry is broken in the z- 
direction. Then the fermion Green's function G(k) becomes l/^Ary + m(k)r 3 ). Note that 
G is a matrix Green's function in both the spinor space and the spin space. The "mass" 
m(k) is related to the sublattice magnetization M by 

d 3 q 



M 



(27T) 



f tr 7iT [G(g)r 3 ]. 



(59) 



Self-consistent equation for m(k) can be obtained from the Schwinger- Dyson equation. Ex- 
pressed in terms of matrix (both in r and 7 space) Green's function and (matrix) self energy 
(pictorially represented by Fig.Ek), the S-D equation is 



E(Jfc) = -m(k)r 3 



d 3 q ~ 



(27T) 



TJk, q)G(k - q^D^q) 



(60) 



where f M , D^ v are fully renormalized vertex and gauge propagator. 

Within the approximation of replacing T M by 7 M and D^ v by D^ u , we have 

d 3 k 'y fl m(k)'j u 



m(p) 



-D^ v (p - k). 



(61) 



(2n) 3 k 2 + m 2 (k) 

Still this is a nonperturbative theory; it is easy to check that a finite order perturbation the- 
ory cannot generate a mass term in the fermion Green's function. Diagrams that contribute 
are shown in Fig.|]b. 



a) 



b) 




FIG. 4. Schematic representation of the Schwinger-Dyson equation. In part a) the solid line 
with shaded blob is the self-consistent Green's function of the fermions G, the thick wiggly line 
is the dressed gauge field (which incorporates the changes in the vacuum polarization due to 
changes in fermion Green's function), and the shaded triangle is the dressed vertex. Part b) is a 
representation of the contribution to the self energy. 
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The context in which the self-consistent equation arises is similar to the SDW (spin 



density wave) problem and the superconductivity ||57|| . However, in our case, the mass 
m(p) is dependent upon 3-momentum. Eq.|61] has been already analyzed by Appelquist and 
coworkers in a different context ("chiral symmetry breaking") [[HJ]. Some of the steps are 
sketched below. 



From the result 



l»lv U M „ - I = 7m7m - = 3-1 = 2 (62) 

we have, after some angular integrals, 

Note that the lattice origin of our theory sets the UV cutoff scale A, while in the theory of 
Appelquist et at which retains the kinetic term ~ -^F 2 U , the coupling constant e 2 sets the 
scale (QED3 is a superrenormalizable theory). The integral equation (Eq.^3]) is equivalent 
to the differential equation 

d_ f 2 dm(p) \ = 8 p 2 m(p) 

dp \ dp J n 2 Np 2 + m 2 (p) 



with boundary conditions 



and 



Am'(A) +m(A) = (65) 



< m(0) < oo. (66) 

It turns out that this nonlinear differential equation has a nontrivial solution only for N < 
N c = 32/7T 2 [KH]. For the physical case of SU(2) antiferromagnet, N equals 2; therefore, the 
dynamical mass generation occurs, and the Neel-vector rotation symmetry is spontaneously 
broken. Thus, provided that we include the gauge fluctuations, we have a Neel order. 

The foregoing argument, however, should be taken with a grain of salt. In principle, the 
gauge propagator that enters the Schwinger-Dyson equation must be a fully dressed one, and 
so should be the vertex. In the symmetry-broken phase, the gauge propagator is different 
from the symmetric phase, as the polarization function is different. At the crudest level, 
if we assume that the fermions acquire a constant mass m(p) = m, then the polarization 
function would be (See Appendix A). 

iWg) ~ -q 2 (<V - , (67) 

m \ q z ) 

which means that a kinetic term ~ ^F 2 U is generated. This results in the Coulomb potential 
for the fermions which in 2+ ID is 

V(x) = -m J d 2 c^- = mm r-^j-j , (R ~ 1/m). (68) 
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Since the potential increases at large distances, it is a confining potential (Sec. V.A), and a 
different physical picture for the fermions might be expected in that case. 

Arguments can be made that these considerations do not affect seriously the conclusion 
as far as the issue of determining whether dynamical mass generation occurs and (in the case 
of occurrence) the value of N c |59|. (The point is that very near N c , the polarization of the 
fermions in the symmetry-broken state must be pretty close to that of the massless fermions. 
On the other hand, above treatment is too crude to study the behaviors of quantities like 
m(p). Lattice gauge theory simulation [^0| does find that the symmetry breaking occurs, 
and iV c fa 3.5, which is close to above analytical results. 

Having seen the symmetry breaking, we can identify the elementary excitations — the 
Goldstone bosons. We all know that the Goldstone bosons in an ordered antiferromagnet are 
spin waves. In the fermion picture, the Goldstone bosons are a collective mode, and appear 
as a pole in the two-particle Green's function, the scattering amplitude in the appropriate 
channel, and in the related vertex [oT|. Here we show this by considering the SU(2) vertex. 
The Bethe-Salpeter equation for the vertex is given by 



A l (p; q) = 1 • t 1 - J-^- 3 ^G(k)A l (k; q)G(k + q) lv D^(p - k), 



(2k) 



(69) 



which is represented by the ladder diagrams of Fig.[5] (In Eq.|69|, the 1 in the first term on 
the RHS is the unit matrix in the spinor space.) 





+ 




FIG. 5. Bethe-Salpeter equation for the isovector vertex in the Q = (tt, it) channel. The fermion 
lines willed a shaded blob represent the renormalized (self-consistent) Green's function. 

If there is a Goldstone boson, then A l (p;q) has a pole at q 2 = 0, in which case the 
homogeneous equation 

A'(p;0) = - 



d 3 k 



7 M G(fc)A i (fc;0)G(fc) 7 ^(p-fc) 



has a nontrivial solution. It is easy to see that 

A l (p;0) = [m(p)r 3 } T l ] 



(70) 



(71) 



with m(p) given by Eq.^T] is the solution. Therefore, if there is a dynamical mass generation 
(m(p) ^ 0), then we have two Goldstone bosons (SU(2) symmetry breaking) : t 1 = t\t 2 
(or r + , r~) in Eq.[7l] gives a nonvanishing commutator. The Lorentz invariance restricts the 
mesons to have a linear dispersion q$ = q 2 (Minkowski space), which is indeed the case with 
antiferromagnetic spin waves. 
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IV. TWO-DIMENSIONAL UNDERDOPED CUPRATES 



A. antiferromagnetic correlations 

The foregoing may appear to be an unnecessarily roundabout way of looking at 2d 
quantum antiferromagnet, yet we believe it is not without certain value. The qualitative 
idea that a strong attraction between spinons and antispinons via gauge field can result in 
the formation of a vector condensate with Q=(tt,tt) (the antiferromagnetic channel) may 
shed some lights on the underdoped cuprates. In the underdoped cuprates, an effective 
theory based on the sFlux ansatz of the SU(2) mean field theory [f)4,l4|,17 consists of 2 



flavors of massless Dirac fermions and a U(l) gauge field just like above, but now the gauge 
field is also coupled to the bosons (holons). In other words, schematically, 

L = ifcHpidn - ia^ipa - ia^jj^ + L B . (72) 

This additional coupling to the bosons will weaken the gauge field in the sense that it will 
screen the time component of the gauge field (lim^o 7^ 0). The gauge propagator then 
would take the form (in the Coulomb gauge) 

H ~ -^o^6 v0 — (73) 




q 2 J x j + v<r q 

where i,j = 1,2 are spatial indices and x is the concentration of doped holes. In the simplest 
approximation, we ignore the massive part, and consider only the spatial component 

DM = (*« - q ~f) , (74) 

and examine the self-consistent equation for m(p) (Eq. |6l|) . 
Since 

l^lvD^ = "fai — = 2-1 = 1, (75) 



we have 



m W = ^V t dk ^rrl \,A k + p-\k- P \). (76) 



This is identical to the Eq.pli except for a factor of 2 difference in the prefactor. This factor 
2 reduction can be understood from that fact that the gauge field in 2d has one transverse 
mode and one longitudinal mode the latter of which becomes massive. From the analysis 
following Eq.|61|, we know immediately that there will be a symmetry breaking only for 
N < N' c = N c /2 = IG/n 2 . Thus, for the physical case of iV = 2, the spontaneous symmetry 
breaking would not occur! — This is what was hoped for our mean field theory. 

The attraction in the Q = (tt, it) channel mediated by the gauge field, although not strong 
enough to generate a condensate, will nevertheless have a strong effect on the spectrum of 
antiferromagnetic excitation. The fluctuation of the order parameter associated with the 
transition (staggered moment) can be examined by looking at the staggered-channel spin 
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correlation function in the ladder approximation, similar to the more familiar problems like 
the superconducting fluctuations or the ferromagnetic spin fluctuations (See Fig.^j) . 




a) b) c) 

FIG. 6. Ladder diagrams, a) Our case: staggered spin correlation. The wiggly lines are 
interactions mediated by the gauge field. Structurally similar examples: b) ferromagnetic spin 
correlation. The dashed lines are short-range repulsive interactions, c) superconducting correlation. 
The dotted lines are some kind of attractive interaction causing pairing. 

In the problem of nearly ferromagnetic Fermi liquids, short range interaction between 
fermions are often modelled in terms of an on-site repulsion U (UN(Ep) is the dimensionless 
coupling constant corresponding to our 1/N). This problem is a lot simpler, as the ladder 
series sums immediately to the RPA form x = Xo/ (1 — Uxo)- The pole of the RPA propagator 
gives the diffusive mode ( "paramagnons" ) associated with a conserved order parameter. This 
mode (more accurately the peak in x"(c<;)) comes down in energy and becomes sharper as 
U approaches U c . 

Unfortunately, in our problem the interaction is retarded and long-ranged, hence the 
diagrams are not easy to evaluate. Nevertheless, on physical grounds, it is quite reasonable 
to expect that the same gauge field which caused the antiferromagnetic instability in the 
absence of holons will try to create a (massive) mode (particle-hole bound state in the 
Q = (jr, tt) channel) in this case. Because the symmetry is unbroken, there is a particle-hole 
continuum, as a result of which a sharp mode cannot exist, but a "broad resonance", i.e. a 
very unstable meson with a Minkowski-space dispersion 

q 2 = rh 2 — irhT, (77) 

is expected. The mass of the mode rh would come down as the transition is approached ( "soft 
mode"). A physical consequence will be that the dynamic susceptibility XQ=( 7r7r )(^) (which 
can be probed by neutron scattering) has a broad peak, with a substantial rearrangement 
of the spectral weight compared to the mean field prediction. This heuristic picture is 
consistent with experiments of Keimer and collaborators [f22|] that find in the normal state 
(and in the superconducting state) of underdoped cuprates a magnetic scattering with a 
broad peak at some frequency scale that comes down in energy as the doping is reduced. 
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FIG. 7. Rearrangement of spectral weight in Xq = ( 7I - 7r )( Ll -') due to gauge fluctuations. 

In principle, we should be able to study this by looking into the Bethe-Salpeter equation 
for spinon-antispinon scattering amplitude or the associated vertex. The analysis, however, 
entails a number of practical and technical difficulties. We can get around the problem 
of the complicated Lorentz-noninvariant gauge propagator for the doped case (Eq.[T3|) by 
convincing ourselves that the we would obtain a similar physics by considering a Lorentz 
invariant one 



D 



fJbV 



8 



(78) 



with N > N c (Reducing the doping would correspond to reducing N). Still, the fact that we 
are investigating a mode that is massive and damped causes complications unseen in the B-S 
equation for the Goldstone bosons of the ordered antiferromagnet. Generally, the scattering 
amplitude and the vertex have many components ("invariant amplitudes"), e.g., 



A z (p; q) = /(p; q)r l + g(p; g)pir l + 



(79) 



The Goldstone modes are massless, hence we can focus on q = in which case /(p; 0) 
decouples from other amplitudes and has the same equation as that of the dynamical mass, 



as we have seen in Sec. III. C. The decoupling doesn't occur for q ^ 0. Appelquist et al. [66 



have considered the scalar component of the Euclidean scattering amplitude in the ladder 
approximation, ignoring the coupling to other components; this corresponds to considering 
f(p; q) only. They claim to find no light meson whose mass comes down to zero as the 
transition is approached, and hence the transition must be a novel one, i.e. neither a first 
order nor a second order transition. However, this conclusion is questionable, as a well- 
defined mode (pole in the scattering amplitude) in the symmetric phase is too stringent 
a requirement for a second order phase transition; more generally what happens in the 
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symmetry-unbroken phase as the second order transition is approached is a build-up of the 
low energy spectral weight for the response related to the order parameter. In any case, the 
analysis of Ref. |)6| is not appropriate for finding an unstable meson which resides in the 



second Riemann sheet of the Minkowski space. 

We now discuss briefly the NMR relaxation rates, which measure the very low energy 
limit of the spin excitation spectrum. As noted earlier, the mean field theories predict 
similar behavior of the Oxygen and Copper site relaxation rates, while experiments show 
marked difference. It would be difficult to calculate finite temperature properties within our 
gauge theory, but we believe that what has been discussed so far throws some light on this 
issue: once the gauge fluctuation is included, there is no reason why the two relaxation rates 
should have similar behaviors, since the effect of the gauge field on the Q=0 and Q=(n, ir) 
response is very different. It is plausible that at high temperature side the Cu \jT\T will 
be increasing with decreasing temperature, as the gauge fluctuation restores certain amount 
of antiferromagnetic correlations (note that above some energy scale, the screening effect 
of bosons won't be too important and the propagator would look like that of the undoped 
cuprates), while at low temperatures (below some scale) it will go down with decreasing 
temperature, since we still have the d-wave gap if the mean field picture is going to survive. 





a) b) c) 

FIG. 8. Possible contribution from boson- fermion coupling, a) The inner loop (dashed line) 
and the outer loop (solid line) denote bosons and fermions, respectively, and the dotted lines 
denote boson- fermion interaction b) "electron" propagator c) schematic representation of the same 
diagram 

Finally, we note that in our discussion of antiferromagnetic excitations in underdoped 
cuprates we have not considered certain class of diagrams that may be important. An exam- 
ple, shown in Fig.|8|, involves a boson loop and a fermion loop exchanging a residual attractive 
interaction. The set of a boson and a fermion line connected by ladder type interaction can 
be viewed as a electronlike propagator (Bethe-Salpeter equation for electron propagator). 
This type of contribution may be related to pieces of the Fermi surface near (±7r/2, ±7r/2) 
fl~4| , and might provide some clue to the origin of low-frequency incommensurate features 
observed in the neutron scattering. 



B. thermodynamic properties 

In the underdoped cuprates the gauge field interaction also affects the uniform part 
of the spin response (the thermodynamic properties) though the effects are subtler. More 
specifically, the coupling of the gauge field to nonrelativistic bosons results in the logarithmic 
renormalization of the velocity of the fermions, which has the effect of enhancing the specific 
heat and the uniform susceptibility ||64j| . This seems to be in accordance with experiments 
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on underdoped cuprates. Although one might alternatively view the enhancement features 
in the normal state as having to do with the "Fermi surface segments" , these are not really 
quasiparticle states in the strict sense (z=0) and there are likely to be theoretical complexity 
and possibility of overcounting. It seems a lot simpler to view that the fermions account 
for most of the entropy and spin response. A curious feature |67| of our theory is that the 
Wilson ratio W = C/Txu has a value quite close to that of the "quantum critical" phase of 

58]: W in our theory is 0.128, while the O(N) nonlinear sigma 
0.124 at zeroth order, and W = 0.116 with the inclusion of 1/N correction. 



the nonlinear sigma model 
model gives W 
Could this be a coincidence? 



V. CONCLUDING REMARKS 

We now conclude with a discussion of some difficult matters, hoping not so much as to 
resolve them but to content ourselves with placing them in perspective. 



A. confinement, spinons, and all that 

The term "spinon" has been used rather carelessly in this paper. In the strict sense, 
it refers to well-defined neutral elementary excitations carrying spin 1/2. A well-known 
(and perhaps the only widely-accepted) example is the solitonic (topological) objects in spin 
chains that are usually understood within the bosonization framework p4| , |26"|| . 

The "spinons" in this paper refer most of the times to spin 1/2 fermions strongly cou- 
pled to fluctuating gauge fields — the fermions that arise from the RVB mean field theories. 
These spinons are far from being well-established. In Id, the fact that they are different 
from topological spinons ( "true spinons" ) has been emphasized by Mudry and Fradkin pS 



Nevertheless, we have seen that in the Id Heisenberg model, the mean field fermions give a 
reasonable description of the specific heat and uniform susceptibility the Zeeman split- 
ting of the spin up and down fermions in a magnetic field and the associated incommensurate 
magnetic excitations |7(J can be easily understood in terms of the mean field fermions. Fur- 



thermore, we showed in Section II that the perturbative treatment of the gauge fluctuations 
systematically improves the antiferromagnetic correlation. Thus the mean field fermions are 
not as bad a starting point as some would regard. 

In 2d and 3d, pre-high T c empiricism and experiences provide strong resistance to the 
notion of spinons as derived from RVB mean field theory. After all, in most metals (2 and 
3d), the basic elementary excitation is quasiparticles that carry both spin 1/2 and charge e, 
and in most insulating antiferromagnets or in spin liquids (like spin ladders), the elementary 
excitations are spin 1 objects (spin waves or magnons). In other words, spinons and holons 
do not exist on their own, but are "confined" in spinon-antispinon composite objects (spin 
waves: /*/) or in spinon-holon composite objects (quasiparticles: f*b). 

The sense in which this "confinement" is discussed is similar to the quark confinement 
in particle physics — the absence outside the nucleus of the quarks that make up hadrons 
(baryons and mesons). In fact, the strong interaction physics appears to share quite a 
few parallels with our problem |7I |. Quantum chromodynamics (QCD), which is widely 



believed to be the correct theory of strong interactions, is a gauge theory whose low energy 
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physics is as poorly understood as the high T c cuprates. At a more substantial level, the 
phenomenology of strong interactions indicates that the (approximate) chiral symmetry is 
spontaneously broken, giving rise to mesons (such as pions) which are Goldstone bosons, 
in analogy with spin waves in a Neel ordered system. Just as the nonlinear sigma model 



gives an excellent description of the low energy spin dynamics of the undoped cuprates |72 



in particle physics it has been well known that similar effective lagrangians (sigma models, 
chiral perturbation theory J75 |, etc.) give a very good description of the hadronic physics. 
However, attempts to "derive" the parameters of the effective theories, such as the pion 
decay constant, from first principles have not been entirely successful. The basic difficulty 
is that the same asymptotic freedom that led to the confirmation of the quark picture at 
the high energy side causes grave difncuties in analyzing the low energy physics. At present, 
no consensus exists as regards the "mechanism" of chiral symmetry breaking in QCD, but 
at least it seems clear that the problem is intimately connected to the general issue of 
confinement. 

It is the basic underlying idea of this paper that the confinement does not occur (at least) 
for the normal state of the superconducting cuprates, or more loosely, spinons and holons 
are "meaningful" objects. On the other hand, some sort of confinement is relevant to the 
discussion of the Neel ordered state in undoped cuprates, since we know that there the low 
energy excitations are not the spinons but spin waves. In considering the confinement, we 
are (helped but also) burdened by previous studies in particle physics regarding the issue. 
The confinement motivated by the strong interaction phenomenology — no color nonsinglet 
particles exist in the physical spectrum — is a strong statement, and it is not clear to what 
extent the confinement in our case should match the quark confinement. 

For example, we can ask the following questions in 2d undoped systems. 1) On the low 
energy side, are the spin waves the only massless degrees of freedom, or could there be an 
additional mode hiding? 2) On the high energy side, do the spinon states exist, or have they 
disappeared completely? In other words, can we somehow "see" spinons at high energies? 

To address these questions, we need to discuss one possibly important aspect of our 
gauge theory, the compactness, that has been ignored so far. The models that we have 
examined in this paper originate from the lattice, hence are compact gauge theories. The 
usual assumption is that the compact theory can be replaced by a more amenable noncom- 
pact theory in the continuum, but this is not always justified. A representative and scary 
example is the pure gauge theory in 2+1D. Polyakov ]74]] has shown that the compact pure 
gauge theory (2+1D electrodynamics) differs from the noncompact one due to instantons - 
the topologically nontrivial, extended classical solutions of Euclidean gauge field equations, 
that can be viewed as tunneling events between topologically inequivalent vacua. Instantons 
cause the Wilson loop to follow the area law (exp(/ dx^a^)) ~ exp(— Area), which means 
the presence of a linear potential V(x. a ,x. b ) ~ |x a — x 6 | between static sources, a sign of 
(strong) confinement. 

When matter fields are present as in our case and in QCD, the situation is not so simple. 
In QCD, for example, despite theoretical attempts |77|] the relevance of instantons to quark 
confinement remains uncertain. Generally, the fluctuations of matter fields (especially the 
massless fields) are adverse to instantons. One specific scenario by which instantons are 
suppressed is fermion zero modes. For example, in the massless Schwinger model, the fermion 
determinant det[(<9 M — 20^)7^] in a gauge configuration with an instanton vanishes, so only 
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the topologically trivial sector contributes to the functional intergral (partition function). 
Similar mechanism turned out to account for the famous U(l) problem [78j in QCD. In both 
cases (1+1D and 3+1D), the zero mode is connected to (axial) anomalies, whose analogue 
in odd space-time dimensions (2+1D) is nonexistent. 

In the context of our problem (compact 2+1D gauge theory with Dirac fermions), 
Marston [21| studied the possibility of fermion zero modes, and concluded that zero modes 
do not exist, suggesting a possible relevance of instantons. Marston has also calculated the 
action of an instanton, and found that it is logarithmically divergent with a prefactor pro- 
portional to the number of flavors (5* oc NlnR). Kosterlitz-Thouless type argument then 
indicates that below critical M c which turns out to be 0.9, instantons may proliferate, while 
for N > Af c which includes the physical case (N = 2), the instantons are suppressed. Ioffe & 
Larkin [^6[ have also calculated the action of an instanton, considering only the bilinear part 
of the gauge field action obtained by integrating out the fermions, and found logarithmic 
divergence, but Af c in this case turns out to be 24. It has been noted that this treatment 



may be too simple, and the question of instantons for physical N remains unclear [79 



The suppression of instantons may no longer be the case if the fermion masses are dy- 
namically generated by spontaneous symmetry breaking, which is indeed our situation. We 
now have the possibility of symmetry-breaking-induced instantons. These instantons here 
cannot induce the dimerization (valence bond solid order) ||8Q|| , unlike those of the bosonic 
spinon theory [jSlH . Probably the most notable consequence of the instantons is that the gauge 
field will become massive J74]|, hence the spin waves will be the only low energy excitations 
of the ordered antiferromagnet. 

On the other hand, in the (simpler) scenario without the instantons (Sec.III.C), there is 
a massless gauge field in addition to the spin waves in the ordered antiferromagnet. This is 
rather bothersome, since it is commonly believed that the spin waves deplete the low energy 
excitations in 2d quantum antiferromagnet. Practically, it may be difficult to determine the 
absence or the existence of the massless gauge field, as it does not couple to external probes 
in a simple way. The gauge field may have a ~ T 2 contribution to the specific heat like the 
spin waves, but the modification of the prefactor due to gauge field might be considered a 
spin- wave renormalization effect. 

The pictures with or without instantons may also have different implications for the 
spinon confinement. If the instantons are relevant, the Polykov-type nonperturbative mech- 
anism might lead to the final picture in which the spinons have vanished completely, and 
traces of them might not be present even at high energies. In the picture without instantons, 
spinons may be still confined, due to the long range nature of the gauge field. A situation 
in which this works out |7^] is the 1+1D CP N ~ l model (L = {(d^ — ia^)z\ 2 + m 2 z Jf z), 
in which integrating out the matter field (z-field) generates a kinetic term (self energy) 
for the gauge field ^F 2 U that gives rise to the linearly confining Coulomb potential (i.e. 
V(xi) = — / dqi exp(iqiXi) / qf ~ \xi\) between the z-charges. In our 2+1D problem 
(Sec.III.C), at a crude level, the dynamical mass generation will have a feedback effect 
on the gauge field, softening the gauge propagator from l/^/q 2 to m/q 2 , making it look like 
a genuine electromagnetic field. In 2+1D, the Coulomb potential associated with electro- 
magnetic field is weakly (marginally) confining (l^(x) ~ log(|x|/_R)) |76]] , and for distances 
{energies} shorter {higher} than some scale R ~ 1/m {m}, the spinon might be "visible". 
Within the framework of Sec.III.C, to see if there are well-defined spinons, we can study 
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whether the Minkowski-space Green's function of the spinons has a pole at some mass scale. 
If the dynamical symmetry breaking occurred in such a way that the spinons acquire a 
constant mass as in the NJL model [p2]| , then the spinons should exist at high energies, 



and we might be able to "see" them |82[ by a method analogous to the Id example |70 
If the symmetry is broken by (long-range, retarded) gauge field interactions as in our case, 
the problem is quite complicated; we need to continue the mass function m(p) from the 
Euclidean space to the Minkowski space (p 2 — > —p 2 ), and examine whether p 2 — m(—p 2 ) = 
has a solution near the real axis. The results from the literature |59(| indicate that when the 



"feedback" effect of the spinon mass generation on the vacuum polarization is taken into 
account (in which case the gauge field dynamics looks like F 2 U ), no poles are found near the 
real axis. 

At present, it is not clear which picture (i.e. with or without instantons) of confinement 
is realized for the undoped cuprates. In the normal state of the underdoped cuprates, in- 
stantons are probably not relevant, as we have massless fermions and there are additional 
massless fluctuations due to the introduction of holes. In any case, the unusual phenomenol- 
ogy of high T c cuprates points to that the spinons and holons are deconfined, as emphasized 
earlier. It is beyond the scope of this paper to discuss the issue of possible confinement in 
the superconducting state. 



B. loose ends 



In this paper we have examined the magnetism of undoped and underdoped cuprates from 
the point of view of neutral fermions with spin 1/2 |83fl . Admittedly, the theory as it stands 
is far from rigorous. The philosophy has been to analyze possibly the simplest effective 
field theory of massless Dirac fermions, bosons and gauge fields, motivated by the sFlux 
phase that appears as a saddle point solution of the SU(2) mean field theory. In reality, the 
situation is a lot more complicated. The mean field fermion spectrum is anisotropic: e(k) = 

\fvpkf + v%k%, and the velocities vp,v 2 have some doping and temperature dependences. 
These dependences, however, do not account for the puzzling properties of the cuprates that 
we have discussed, and the photoemission does indicate that the quasiparticle gap remains 
large in the normal state of the underdoped cuprates (i. e. superconducting transition is not a 
gap-closing transition) and the gap is only weakly doping dependent. Therefore, the effective 
theory may be quite sensible for studying qualitative features not captured by mean field 
theories. We have made a number of approximations to treat the complicated dynamics of 
the gauge fields, but we hope it's not too optimistic to view that the qualitative conclusions 
are correct. In any case, the idea that the holon coupling to the gauge field prevents 
spontaneous symmetry breaking in the fermion system is an attractive one, and is very much 
in the spirit of the empirical fact that moving holes quickly destroy antiferromagetic order. 
Unfortunately, even with simplifications, the calculations quickly become rather intractable. 

Eventually the "spin gap" has to close up as we go to the optimally doped regime, which 
in the mean field theory is modelled by the uRVB saddle point |[53|| . The details of this 
crossover is certainly beyond our hopes. This will severely limit us in considering some very 
interesting issues, like the relation between the neutron scattering peaks in the underdoped 
cuprates and the sharp 41meV peak in the superconducting phase of the optimally doped 
YBCO7. Anyway, the present theory does not say much about the spin excitations in the 
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superconducting state of the underdoped cuprates, since the dynamics of gauge field in the 
superconducting state is different from that considered here. Another drawback is that in 
our framework it is not easy to consider possible incommensurate features in spin excitations. 

Granting these limitations, we still feel that the picture of spin excitations in underdoped 
cuprates in terms of deconfined fermionic spinons is reasonable and perhaps more natural 
than other descriptions like those based on fluctuating staggered moments. Features like the 
linear-in-T behavior of uniform susceptiblity and the specific heat coefficient in underdoped 
cuprates might be also explanable in terms of the "quantum critical" regime of the nonlinear 
sigma model, but it is not clear in that approach how to account for the strange behavior of 
the Copper 1/T{T in the same temperature range. Attempts to explain the Copper 1/71 T in 
terms of the "quantum disordered" regime then has to explain why activated behaviors are 
not seen in quantities like uniform susceptibility. Again, this seems to point to the difficulty 
of achieving a theoretical description of a system that involves a mysterious combination of 
gaplike (short range) correlations and critical correlations. 
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In this appendix, fermion vacuum-polarizations in 1+1D and 2+1D are worked out using 
"dimensional regularization" . 

For massless Dirac fermions, the polarization function U^ u is given by 
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APPENDIX A: FERMION POLARIZATION 




= N 



f dx 



d D k' tr[(fc' + (1 - x)q) 11 ^k' - xq)^ v ) 
(2n) D {k 12 + q 2 x(l - x)) 2 



(Al) 



where N is the number of flavors of fermions. Using the trace identity 



3 puffier ) ) 



(A2) 



The numerator of the integrand is found as 



2k'X - 2x(l - x)(q^ - 5^q 2 ) - 5^[k' 2 + q 2 x(l - x)]. 



(A3) 



Substituting, we get 
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2k ^k v 



(27 

Sftv 2x(l - x)(5^ u q 2 - q^q, 



k 2 + q 2 x(l — x) (k 2 + q 2 x(l — x)) 2 
The first two terms cancel, and we have 



(k 2 + q 2 x(l — x)) 2 

(A4) 



U, u (q) = 2AT(trl) r(2 f ; ( 2 V <W ~ %<1») jf dxx(l - x)(q 2 x(l - x)) ' 2 ' 2 



(Au) D / 2 

^yff^-tf) D = 2 + l 
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m^if 5 u -!Bf) d = i + i. 



(A5) 



For the situations discussed in the main text, we have trl=4 for the 2+1D (Sees. Ill & IV) 
as we have 4 component fermions, while trl=2 for the 1+1D (Sec. II). 

If the fermions are massive, i.e., L = ^(d^ — ia^^ip + mtpip, similar calculation shows 
that the polarization function is given by 

U^q) = 2N(tTl) T ^~^l 2 2 \ q\ u - q,q u ) jf* dx(l - x)x(m 2 + q 2 x{\ - x)) ' 2 ' 2 

= %» - W.) [ 7 - ? ^ ?T? tanh'( i ^) ) [D = 1 + 1], 

(trlW, r. . ( m q 2 — 4m 2 q \ r 

= - wo + ^^^-(^=1=5) J [D =2 + i]. 

For small g (i.e. q 2 « m 2 )), we would have 

D = 2 + l:U llu (q)^^^(q\ u -q ll q u )(^-- ''~ ^ 



47r y6m 60m 3 

D = 1 + 1 : rV(g) » ^^(q\ u - q^ ( 1 



2tt v " ^ \6m 2 30m 4 ' / 

Note that the massive fermions give rise to a gauge field dynamics of the form L ~ alia ~ 
{F^ u ) 2 like the real electromagnetic field. 

APPENDIX B: JACOBIAN FOR CHIRAL ROTATION 

In this appendix, we derive the jacobian exp(J) = for chiral rotation, along 

the line of Stone |35|. Basically we require the functional integral Z — J DipDip exp(— L) 
be invariant under the chiral transform. For an infinitesimal chiral transform if) — > if)' = 
exp(— 75# a ) , 0) the measure changes as 

D${x)Dij){x) ^D$(x)Di>(x) (l-J d 2 z|^0 o (x)^ , (Bl) 
while the lagrangian changes as 
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L(VO -> L(VO - {dMwrt = L(ip) + z^a M j 5/1 . (B2) 
Setting the functional derivative 

f = " ia *W (B3) 

to zero gives 

w. = -s^ F - = + 7 fl ! s " (B4) 

where we have used the axial anomaly condition <9 M j 5M = —^e^ u F^ v with I 7 ^ = d^ — d^a^, 
a>n = ^iyd v 9 a + d^6 b . Eq.[B3] then implies 

r N 

J=~ d 2 x—(d,9 a ) 2 . (B5) 

J Z7T 
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